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EXECUTIVE  SUMMARY 


Sound  absorption  is  a  major  consideration  in  SONAR  systems. 
Absorption  in  sea  water  involves  pressure-dependent  chemical  equilibria  with 
relatively  slow  relaxation  rates  Understanding  the  mechanism  is  essential 
to  account  for  effects  of  variable  environmental  factors. 

Boric  acid  was  identified  as  the  cause  of  low-frequency  absorption 
by  means  of  a  temper ature- jump  (T-jump)  method  which  measures  rate  of  change 
of  pH.  Concentration  studies  showed  that  the  relaxation  frequency  never 
exceeds  4  kHz  in  either  fresh  water  or  synthetic  sea  water.  Results  in  both 
cases  are  consistent  with  a  rapid  hydroxyl  association  followed  by  the  rate¬ 
controlling  monomolecular  step.  However,  theoretical  absorption  predictions 
are  much  too  low  in  either  case. 

Absorption  measurements  in  synthetic  sea  water  by  the  resonator 
method  now  show  that  the  relaxation  frequency  can  not  only  exceed  4  kHz,  but 
also  depends  on  the  concentration  of  carbonic  acid  as  well  as  boric  acid.  In 
addition,  absorption  depends  on  the  concentration  product.  These  results  are 
consistent  with  a  slow  single-step  exchange  reaction  which  does  not  involve 
change  in  pH.  Coupling  to  the  fast  CaC03  ion-pair  equilibrium  is  also 
important  because  absorption  is  too  small  if  calcium  is  omitted. 

The  discrepancy  between  resonator  and  T-jump  measurements  is 
resolved  if  both  steps  of  the  hydroxyl  association  are  slow.  In  sea  water, 
the  coupled  system  then  has  two  slow  relaxations,  the  association  controlling 
pH  and  the  exchange  controlling  absorption. 

It  is  concluded  that  the  dominant  low-frequency  absorption  mechanism 
in  sea  water  is  the  exchange  reaction  and  that  the  laboratory  measurements  of 
synthetic  sea  water  are  completely  consistent  with  the  empirical  formula 
derived  from  sea  measurements. 
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ABSTRACT 


Low-frequency  sound  absorption  in  sea  water  is  dominated  by  a 
1  kHz  relaxation  involving  boric  acid.  Resonator  measurements  in  synthetic 
sea  water  show  that  the  relaxation  mechanism  can  be  modeled  as  an  exchange 
between  the  borate  and  carbonate  equilibria  in  which  coupling  with  calcium 
ion-pair  equilibria  provides  most  of  the  volume  change.  Consistency  with 
the  empirical  pH-dependent  absorption  formula  derived  from  ocean  measurements 
is  also  demonstrated. 


BACKGROUND 


The  1935  Naval  Research  Laboratory  (NRL)  report  by  Stephenson1 
appears  to  be  the  first  to  recognize  that  sound  absorption  in  the  ocean  is 
much  greater  than  in  fresh  water.  However,  it  was  not  until  after  World 
War  II,  when  the  anomaly  had  been  verified  by  extensive  sea  measurements, 
that  efforts  were  undertaken  to  determine  the  cause.  Laboratory  measure¬ 
ments  showed  that  fresh  water  absorption  follows  the  classical  frequency- 
squared  dependence.  The  magnitude  was  found  to  be  somewhat  greater  than 
predicted  from  shear  viscosity,  but  this  was  eventually  explained  by  bulk 
viscosity.  The  sea-water  anomaly  was  more  than  an  order  of  magnitude 
greater . 


In  1948,  Liebermann2  reported  field  measurements  comparing  fresh 
water  and  sea  water  absorption  in  the  common  frequency  range  of  10  -  1000  kHz 
Earlier  sea  data  were  concentrated  in  the  SONAR  range  below  100  kHz,  while 
fresh  water  had  been  measured  at  much  higher  frequencies.  Measurements  in 
a  reservoir  showed  agreement  with  the  earlier  fresh  water  results.  Sea 
measurements,  on  the  other  hand,  agreed  with  earlier  salt  water  values  in  the 
lower  range,  but  decayed  to  fresh  water  values  at  high  frequencies.  A  relax¬ 
ation  process  was  indicated. 

Relaxational  absorption  can  occur  in  chemical  solutions  as  the 
result  of  pressure-dependent  equilibria  that  affect  the  complex  compress¬ 
ibility  of  the  medium.  The  rate  of  equilibration  determines  the  relaxation 
frequency  which,  in  this  case,  was  found  to  be  150  kHz.  Liebermann  proposed 
the  ionization  equilibrium  of  NaCl  as  a  possible  cause. 
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In  1949  Leonard,  et  al^  reported  initial  results  of  a  laboratory 
investigation  of  synthetic  sea  water  (SSW)  by  the  resonator  method.  The 
resonator  a  12  liter  spherical  flask,  was  excited  by  acoustic  pulses  at 
selected  resonant  frequencies  in  the  range  of  25  -  500  kHz.  The  sphere  was 
suspended  in  a  vacuum  chamber  in  order  to  minimize  residual  losses  and 
improve  sensitivity.  Relaxational  absorption  spectra  derived  from  measure¬ 
ments  of  decay  rate  showed  MgS04  to  be  the  cause  of  the  150  kHz  relaxation. 

A  more  detailed  paper  by  Wilson  and  Leonard*  appeared  in  1954. 

In  1953,  Kurtze  and  Tamm®  also  confirmed  Leonard's  findings  and,  in 
addition,  investigated  the  interactive  effects  of  NaCl.  A  number  of  possible 
mechanisms  were  proposed.  However,  it  was  not  until  1962,  in  a  definitive 
analysis  of  Eigen  and  Tamm,®  that  the  details  of  the  multi-step  MgSO^  mech¬ 
anism  were  delineated. 

In  1962,  Schulkin  and  Marsh?  introduced  a  formula  for  sound 
absorption  in  sea  water  based  on  the  MgSC>4  results  and  ocean  measurements  in 
the  frequency  range  of  2  -  25  kHz.  When  the  focus  of  attention  in  underwater 
acoustics  shifted  to  lower  frequencies,  the  validity  of  extrapolating  the 
formula  was  subject  to  question. 

In  1965,  Thorp,®  in  a  summary  of  low-frequency  attenuation  measure¬ 
ments,  found  values  consistently  an  order  of  magnitude  greater  than  predicted 
(figure  1-1).  He  also  showed  that  the  excess  could  be  modeled  as  a  1  kHz 
relaxation.®  Similar  measurements  in  the  Mediterranean  by  Skretting  and 
Leroy  showed  a  1.7  kHz  relaxation. 10  A  chemical  relaxation  similar  to  MgSC>4 
seemed  to  be  the  most  promising  explanation. 

To  investigate  the  proposed  second  chemical  relaxation.  Browning, 
et  al11  carried  out  a  series  of  sound  channel  propagation  measurements  with 
different  ambient  temperatures.  A  systematic  increase  in  relaxation  frequency 
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Figure  1-1.  Thorp's  Compilation  of  Low-Frequency  Attenuation  Data 
Compared  with  the  Marsh-Schulkin  Formula  (Second  term  in  the 
equation  for  solid  curve  is  the  1  kHz  relaxation.) 
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with  temperature  would  be  expected.  Fresh  water  experiments  were  included 
for  comparison  purposes.  Reports  and  publications  resulting  from  this 
program  are  compiled  in  reference  12. 

Results  from  Lake  Superior18  showed  low-frequency  attenuation 
comparable  to  Thorp's  anomaly  suggesting  that  sea  salts  might  not  be  respon¬ 
sible.  However,  the  Lake  Tanganyika  experiment1*  gave  no  evidence  of  a  low- 
frequency  relaxation  and  the  hypothesis  was  abandoned.  The  Lake  Superior 
anomaly  was  subsequently  ascribed  to  leakage  caused  by  scattering  from  random 
inhomogeneities  in  the  thermocline.  Frequency-independent  attenuation,  as  the 
result  of  diffusion,  is  apparently  common  to  all  sound-channel  experiments. 
However,  the  loss  is  usually  small  and  detectable  only  at  very  low  frequencies 
because  of  masking  by  absorption. 18 

Lake  Tanganyika  results  a^so  showed  excess  high-frequency  loss 
similar  to  the  MgS04  relaxation  in  sea  water.  This  lake  is  unusually  high 
in  carbonate,  but  magnesium  and  sulfate  concentrations  are  much  too  low  for 
MgS04  to  be  responsible.  Investigation  of  synthetic  lake  water  by  the 
resonator  method  showed  the  cause  to  be  a  similar  MgC03  relaxation. 16  The 
same  relaxation  also  occurs  in  sea  water,  but  it  is  only  a  minor  absorption 
component  compared  to  the  others.17 

Results  of  a  number  of  sea  experiments  were  then  fitted  by  a  Thorp- 
type  equation  plus  a  frequency-independent  loss  term. 118  Good  correlation 
was  found  between  relaxation  frequency  and  temperature,  indicating  that  a 
chemical  relaxation  was  the  likely  mechanism.  Considerable  spread  in  the 
absorption  coefficients  was  also  noted.  Following  the  discovery  of  boric 
acid  as  the  cause  of  the  relaxation,  correlation  with  pH  also  showed  the 
predicted  behavior.18 

Because  of  size  limitations,  the  resonator  method  is  not  practical 
for  investigating  a  relaxation  in  the  1  kHz  range.  In  1973,  Yeager,  et  al1® 
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employed  the  T-jump  method  to  study  the  mechanism.  The  relaxation  frequency 
was  determined  from  the  rate  of  change  of  pH  following  a  sudden  increase  in 
temperature.  Synthetic  sea  water  measurements  showed  a  relaxation  in  the 
1  kHz  range  only  with  boric  acid  present.  An  acid-base  exchange  between  the 
borate  and  carbonate  equilibria  was  proposed  as  the  likely  mechanism,  and 
possible  interactive  effects  of  ion-pairing  were  also  pointed  out. 

In  1975,  Simmons20  reported  a  T-jump  study  of  relaxation  frequency 
as  a  function  of  B(0H>3  concentration.  Results  for  synthetic  sea  water  and 
NaCl  solution  were  both  fitted  by  two-step  processes.  Absorption  measure¬ 
ments  were  also  made  in  the  frequency  range  of  6  -  1 50  kHz  using  a  100  liter 
resonator.  Combined  results  appeared  to  support  the  conclusion  that  the 
mechanisms  are  similar  in  both  H2O  and  sea  water. 

Since  1979,  an  investigation  of  sound  absorption  in  sea  water  by 
the  resonator  method  has  been  carried  on  by  Mellen,  et  al. 21-23  Relaxations 
involving  MgCOj0  and  MgB(0H)+  have  been  identified  and  modeled  as  two-step 
processes  incorporating  interactive  effects.  However,  attempts  to  model  the 
BfOHjj  relaxation  by  a  similar  process  proved  unsuccessful. 

In  this  report,  analysis  of  the  resonator  measurements  is  presented 
which  supports  the  first  hypothesis,  i.e.,  that  the  dominant  process  in  sea 
water  is  an  exchange  between  the  borate  and  the  carbonate  equilibria  in  which 
calcium  ion-pairing  provides  most  of  the  volume  change.  In  addition,  con¬ 
sistency  with  the  empirical  pH-dependent  absorption  formula  derived  earlier 
from  ocean  measurements  is  demonstrated. 
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Section  2 


ANALYSIS 


2.1  SEA  WATER  CHEMISTRY 


The  simplifying  assumption  is  that  the  concentrations  of  major 
sea  water  components  vary  only  with  salinity  and  remain  in  constant  ratio. 
Concentrations  of  chemical  species  can  depend  additionally  on  temperature, 
pressure  and  pH.  If  the  relaxation  mechanisms  are  known  in  detail  and  specie 
concentrations  and  other  parameters  can  be  calculated,  sound  absorption 
should  be  predictable  for  any  ocean  conditions. 

Specie  concentrations  depend  on  the  properties  of  the  various  equi¬ 
libria  involved.  In  Leonard's  experiments,  the  synthetic  sea  water  constit¬ 
uents  were  Na2SC>4  and  MgCl2»  while  the  relaxing  species  turned  out  to  be 
MgS04.  The  multi-step  relaxation  can  be  represented  for  present  purposes  as 
the  two-step  equilibriums 

Mg2+  +  S°42'  :  (MgSO°)j  ♦  (HgSO®)^  ,  (2-1) 

where  step  2  controls  the  relaxation  rate  in  question. 

The  stoichiometr ic  association  constant  is 


KMgS 


10/M 


(2-2) 


where  |MgS04°J  *  |(M9S04°)lJ  +  J ( M9S04°)lI | •  Brackets  indicate  concentra¬ 
tions.  Units  of  M  are  gram  moles/liter.  Values  are  for  25*C  and  atmospheric 
24 

pressure. 


2-1 


TR-335N 


Common  ions  cause  coupling  between  equilibria.  In  this  example, 
the  most  important  interaction  is  the  fast  reaction: 


Na+  +  SO.2-  *  NaSO.-  , 

4  *  4 


for  which  the  association  constant  is 


NaSO . 


NaS 


Na+  so  2~ 

4 


2/M 


(2-3) 


(2-4) 


Equilibrium  "constants"  may  vary  with  ionic  strength  (Debye-Huckle 
effect).  In  conventional  sea  water  stoichiometry,  ionic  strength  is  kept 
constant  (I~0.7  M)  while  the  ratios  of  the  various  components  are  varied. 
This  is  unnecessarily  restrictive  and  very  inefficient  for  absorption  meas¬ 
urements  because  everything  must  be  changed  for  each  experiment.  It  is  much 
mote  convenient  simply  to  increase  the  concentration  of  the  component  under 
study,  letting  I  change.  Deviations  from  theory  at  low  ionic  strengths  may 
result  from  variations  in  activity. 

Not  all  major  components  of  sea  water  affect  B(OH)3  absorption. 
Sulfate  will  be  omitted  throughout  because  the  added  MgSC>4  absorption  is  an 
unnecesary  complication.  The  essential  constituents  and  concentrations  for  a 
synthetic  sea  water  that  satisfies  the  requirements  of  the  present  investi¬ 
gation  are  listed  in  table  2-1.  Acid-base  equilibria  involving  B(OH)3  and 

NaHC03  will  be  of  concern.  The  B (OH) 3/B (OH) 4  equilibrium  shows  no 

-  2- 

significant  effects  of  variable  activity.  On  the  other  hand,  HC03  /C03 
activity  changes,  but  becomes  stabilized  for  I  >  0.1  M. 
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Table  2-1.  Synthetic  Sea  Water  Recipe 
(No  Sulfate) 


Constituent 

Concentration  (mM) 

NaCl 

400 

MgCl2 

50 

CaCl2 

10 

NaHC03 

2.5 

B(OH)3 

0.43 

The  equilibrium  constants  of  acid-base  systems  can  be  measured 
by  pH  variation  as  a  function  of  concentration  of  the  major  pair-forming 
metallic  ions;  namely,  Na+,  Mg^4-,  and  Ca2+.  For  example,  in  aqueous  boric 
acid,  the  only  significant  species  at  concentrations  less  than  0.1  M  are 
B (OH) 3  and  B(OH)4  .  The  equilibrium  constant  can  be  written 


B(OH)4"j 

[H1 

i  J 

5.9  x  10_10M 


(2-5) 


where  is  the  hydrogen  ion  concentration.  Only  the  borate  ion  forms 
pairs.  Therefore,  the  apparent  equilibrium  constant  in  sea  water  can  be 
wr itten 


K  » 


'1  fl 


B(OH) A  I  + 

4 


NaB(0H)4°]  +  [MgB  (OH)  *  j  +  [caB(0H)4+]j 

B(0H)3] 


(2-6) 


In  terms  of  the  association  constants,  equation  (2-6)  becomes 


h  ■  KB  Sb  I"’2*!  ♦  'c.  ic*21)  •  |2-7' 
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The  most  straightforward  method  of  measuring  values  is  to  use  borax 
at  low  concentrations,  e.g.,  [B]~10“J  M.  At  low  concentration,  Na  ion¬ 
pairing  is  negligible  and  borax  dissociation  in  H20  yields 


Na2B407*10H20 


2Na  +  2  B (OH) ,  +  2  B (OH)  +  3H.0  .  (2-8 

i  4  i 


With  [B(OH)3]  *  [b(OH)4  |  in  Equation  (2-5),  we  have  KB  =  [H+]  or  pH  = 

-log  [H+]  =  pKg  =  9.23.  When  NaCl,  MgCl2  or  CaCl2  are  added,  pH  decreases 
because  of  ion-pairing  and  [H+]  increases  to  compensate  for  the  reduced 
B (OH) 4  concentration.  Kg  being  constant.  In  equation  (2-7),  we  then 

*  j.  * 

have  Kg  =  (H+]  or  pH  =  pKg  <  9.23.  The  charge  conservation  condition  for 
changes  in  [H+]  can  be  neglected  in  the  pH  range  of  interest. 


The  ion-pair  association  constants  can  be  determined  by  stepwise 
increasing  the  concentrations  up  to  the  sea  water  values.  The  respective 
equilibrium  constants  are  calculated  from  the  slope  of  K0*  vs.  concentration. 
The  ion  concentration  can  usually  be  approximated  as  the  total  concentration 
of  the  constituent,  e.g.,  [Na+]  a  [NaCl] . 

For  calculating  the  boric  acid  concentration, 


'B<°H>3]  -  J! 


pH-pKo 

where  x  *  1(r  For  the  borate  ion  concentration. 
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The  following  values  are  indicated  for  sea  water  at  25*C:24 
KNaB  =  2. 1/M,  KMgB  =  18/M  and  KcaB  *  26/M.  In  figure  2-1,  the  boric  acid 
SSW  model  is  compared  to  dilution  measurements  in  CC^-free  sea  water  by 
Lyman.25  Agreement  appears  to  support  the  overall  validity  of  the  SSW  model. 

In  the  CO2  system,  there  are  two  equilibria  to  consider 

1  2 

H2  C°3  t  HC03~  +  H+  Z  co32~  +  2H+  •  (2-11) 


The  corresponding  values  are  pKjc~6  and  P^C'"'10-  Only  the 
second  equilibrium  appears  to  play  any  significant  role  in  sea  water  absorp¬ 
tion.  Therefore,  the  first  equilibrium  can  be  neglected  for  pH  >  1.  The 
apparent  equilibrium  constant  becomes 


2C 


=  K 


2C 


,U\a2C 


K«,2C  ("5 


2+- 


+  K 


Ca2C 


[ca2*]' 


(2-12) 


The  following  values  are  indicated  for  sea  water  at  25*C:24  KNa2c  *  1.6/M, 
KMg2C  =  5 3/M  and  KCa2c  =  125/M.  In  figure  2-2,  the  carbonate  model  is  com¬ 
pared  with  Lyman's  sea  water  dilution  measurements  of  K2c-  Note  that  the 
nonlinear  region  near  zero  (dashed  curve)  may  come  from  variations  in  the 
association  constants  as  well  as  the  K2c*  The  stable  value  of  K20 
(horizontal  dashed  line)  is  1.8  x  10-10  or  pK2c  =  9.73,  while  the  actual 
value  is  4.7  x  10-11,  or  pK2c°  =  10.33.  Note  that  the  slopes  refer  to  the 
latter  and  are  greater  than  the  respective  association  constants  by  a  factor 
of  4.  The  stable  values  are  used  for  calculation  of  specie  concentrations 
as  in  equations  (2-9)  and  (2-10). 


In  addition  to  specie  concentrations,  values  of  the  volume  changes 
associated  with  the  various  equilibria  are  needed.  Conventional  units  for 
AV  are  milliliters/gram  mole  or  liters/kilogram  mole  and  will  be  abbreviated 
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Apparent  values  may  be  calculated  from  measurements  of  the  change  of  apparent 
equilibrium  constants  with  pressure  using  the  relation 

* 


* 

AV 


-3  InK 
3  P 


RT 


(2-13) 


where  T  is  absolute  temperature  and  R  is  the  gas  constant. 


Another  method  is  by  measuring  the  density  of  solutions  by  means 
of  the  picnometer.  Ward  and  Millero^  used  this  method  to  study  boric  acid 
ionization  in  H2O.  Measurements  of  B(0H>3  and  NaB(OH)^  solution  at  25*C  are 
shown  in  figure  2-3.  Constituent  in  the  latter  is  sodium  metaborate,  which  is 

assumed  to  dissociate  as  follows:  NaB02*2H20  ♦  NaB(0H)4°  *  Na+  +  B(0H>4  . 

Note  that  the  apparent  molal  volume  of  B(0H>3  solution  is  39.2/kM  independent 
of  concentration  while  ion-pairing  in  NaB(0H)4  solution  causes  a  progressive 
increase.  At  small  concentration,  association  of  NaB(0H)4°  is  negligible  and 
the  apparent  molal  volume  approaches 

V*  |NaB(OH)4j  =  V°^B(OH)4~j  +  V°^Na+ j  .  (2-14) 

For  V°(Na+)  -  -5.7/kM,  V°  |b(OH)4_|  *  20.7  +  5.7  *  26.4/kM.  The  value  of  AV 
for  boric  acid  ionization  depends  on  which  equilibrium  is  involved. 


B(OH)3  +  H20  l  B(OH)4"  +  H+  , 

or  (2-15) 

B (OH) ,  +  OH”  ♦  B (OH)  ” 

The  volume  change  for  formation  of  the  ion-pair  NaB(OH)4°  is  calculated 
from  curve  fit  to  the  data.  For  K)jag  *  2.1/M,  the  estimated  value  is 
AVNaB  ”  8/kM.  Other  values  estimated  from  MgB(OH)4+  absorption  measurements 
are  AVj^gg  a  12/kM,  and  AV^ag  *  0/kM. 

Corresponding  values  for  the  CO2  system  are^4  AVj4g2c  “  13. 8/kM, 
AVNa2C  =*  12/kM,  and  AVCa2c  =  12.5/kM. 
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Apparent  Molal  Volume  (kM 


0  0.5  1-0 

Concentration  (M) 


Apparent  Molal  Volumes  of  Boric  Acid  and  Sodium  Borate  * 


H^O  Solution  vs.  Concentration 
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RELAXATION  KINETICS 


The  simplest  case,  representing  a  change  in  structure,  is  the  mono- 
molecular  equilibrium27 


A  *  B 


(2-16) 


The  rate  equation  is 


dA  +  >  -  , 

^  +  k  A  -  k  B  =  0. 


(2-17) 


where  k+  and  k”  are  forward  and  backward  rate  constants,  respectively,  and 
primes  indicate  instantaneous  concentrations.  For  small  perturbations,  let 


A  =  A  +  AA 
B'  =  B  +  AB 


(2-18) 


where  A  and  B  are  mean  concentrations.  Then  AA  =  -AB. 


Substituting  equation  (2-18)  in  equation  (2-17)  and  using  the  equi¬ 
librium  constant. 


_  k+  B 
K  =  \C  "  A 


gives 


+  k  AA  =  0 


(2-19) 


(2-20) 


A  solution  of  equation  (2-20)  is 


AA  «  AA^xp  (-t/x) , 


(2-21) 


where  x  is  the  relaxation  time.  Then  1/x  *  k+  +  k"  =  u>r  where  wr  is  the 
angular  relaxation  frequency. 


2-10 


TR-335N 


For  a  sinuisoidal  forcing  function  equation  (2-20)  becomes 

dAA 

t  — +  AA  =  Q  exp  (iwt)  r  (2-22) 

where  m  is  the  angular  frequency  and  Q  is  the  amplitude.  In  this  case  the 
solution  can  be  written 


AA  =  Q  exp  (+iwt)/(1+iwT>  .  (2-23) 

An  acoustic  plane  wave  with  attenuation  can  be  written 

Q  =  Q0  exp  (-iwx/c)  =  Qq  exp  (-(a+i<)x)|  ,  (2-24) 


where  c  is  the  complex  sound  speed.  Taking  the  compressibility  as  the  sum 
of  the  normal  compressibility  0O  =  |p0c02)-1  and  the  relaxational  compres¬ 
sibility  Br,  the  complex  sound  speed  becomes 


c 


Bc  /  (1+iwt)) 


-1/2 


(2-25) 


where  p0  is  density.  Expanding  the  square  root  for  small  values  and 
separating  into  real  and  imaginary  parts,  we  find 


ic 


1  + 


2  2 


a 


ti) 

c 

o 


cat _ 

1  +  M2T2 


(2-26) 


The  dispersion  term  is  small  and  will  be  neglected.  In  some  work,  the 
attenuation  coefficient  in  dimensionless  form  as  aX,  where  X  is  the  acoustic 
wavelength.  The  absorption  spectrum  is  then  written  as 


aX  =  (<*X)max 


2  mi  _ 

1  +  u/2x2 


(2-27) 
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This  approach  may  be  useful  for  comparing  measurements  at  different  tempera¬ 
tures.  However,  in  all  but  simple  cases,  it  leads  to  extremely  cumbersome 
equations.  In  this  work,  the  alternative  dimensional  expression 


a 


“max 


u>2t  2 

1  +  o>2 1  2 


(2-28) 


is  used  because  the  equations  are  much  simpler  and  easier  to  understand. 
This  fact  can  easily  be  verified  by  converting  equations  in  to  (“Mmax 

for  the  coupled  reactions  discussed  later. 


The  chemical  compressibility  can  be  expressed  in  terms  of  other 
thermodynamic  parameters.  Simplifying  approximations  include  neglect  of 
enthalpy  effects  which  are  small  when  water  is  the  solvent.  For  the  case  of 
equation  (2-16),  also  neglecting  activity  effects,  the  chemical  compressi¬ 
bility  is 


e 

r 


(AV)  2 
RT 


(2-29) 


where  *V  is  the  change  in  molal  volume,  i.e.,  VB-VA. 


Finally, 


substituting  equation 


(2-29) 


in  equation  (2-26)  and  setting 


2  + 

(AV)  k  Aooco_  (2-30) 

max  2  RT  Np/unit  distance  . 

Note  that  equation  (2-30)  can  be  written  in  a  number  of  ways  using  the 
rel ations 


k+A  =  k'B  =  -^-k—  (A+B) 
k++k~ 


(2-31) 
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From  equation  (2-36),  it  is  apparent  that  the  first  step  acts  like  a  coupled 
equilibrium  slowing  the  forward  rate  of  the  monomolecular  step  by  an  amount 
depending  on  the  degree  of  coupling,  i.e.,  for  Ki (A+B)  >>  1,  the  effect  is 
negligible. 


Similarly,  for  Omax'  t^ie  sam®  equation  applies  as  for  the  monomolec¬ 
ular  case  equation  (2-30),  except  that  AV  now  also  reflects  coupling.  For 
K-j  (A+B)  >>  1,  the  effect  of  AV^  is  a).so  negligible. 


Similar  effects  occur  for  other  coupled  reactions  involving  common 
ions.  For  example,  consider  the  system  with  two  fast  associations  coupled  to 
the  slower  monomolecular  equilibrium 


A  +  B 


1  2 
+  E  *  C  +  E  t 


D  +  E 


3\\ 


A  +  F 


*4 


(2-37) 


One  of  the  steps  in  the  fast  loop  is  redundant  because  AVi  =  AV3+AV4  and 
K1  =  k3k4*  Neglecting  step  4,  the  parameters  are 


and 


+  Kt (A+B)+K1K3B(A+B+E) 

“r  ”  k2  +  k2  1+K1  (A+B )  +K3  (B+E)  +1^ (A+B+E) 

AVi(  1+K3  (B+E))  -AV3K3E 

AV  =  AV  _  +  - 5— — - - - 

2  1+K1 (A+B) +K3  (B+E) +K1K3B (A+B+E) 


(2-38) 


Another  type  of  association  characteristic  of  ion-pairing  for  B(OH)3 
in  NaCl  solution  is 


1  2 

A  +  B  +  E  :  C  +  E  :  D  +  E  .  (2-39) 

4  ft  3  5+t 

F  t  G 
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When  steps  2  and  3  are  both  slow,  and  steps  1,  4,  and  5  are  fast,  there  is 
only  one  slow  relaxation.  The  relaxation  frequency  is  given  by 

k2”  +  k3"  K5(D+E)  K^A+B)  (k2+  +  k3+  K4(D+E))  .  (2-40) 

“r  (D+E)  +  '  1  +  K1  (A+B)  |1  +  K4(D+eT)— 

The  exchange  or  transfer  equilibrium  will  be  written^ 

A  +  B  t  C  +  D  .  (2-41) 


In  this  case,  the  parameters  are 


and 


uir  =  k+  (A+B)  +k“  (C+D) 

(AV) 2k+ABp0c0 
°max  2RT 


(2-42) 


The  equilibrium  constant  is  given  by  K  *  k+/k  =  CD/AB.  For  acid- 
base  systems,  we  can  also  have  the  individual  ionization  constant,-  e.g., 

Kt  =  C[H+]/A  and  K2  =  B(H+]/D;  then  K  =  K-\/K2’ 


The  coupled  equilibria  arising  from  ion-pairing  will  be  written 


1 

A+B+EtC+D+E  (2-43) 

3  +t  2  t+  4 

B  +  F  t  C  +  G 

where  steps  3  and  4  are  both  fast.  The  association  constants  are  K3  *  F/AE 
and  K4  =  G/DE.  For  this  system,  we  have  only  one  slow  relaxation  and  the 
parameters  are 
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u>r  =  k*  A  +  k~  D  +  k*  (F+B)  +  k”  (C+G) 

1  k -j  ~  -  k  2~  )c(  1+K3  (A+E)  +K$d)  +  (  k  1 +  -  k2+/  b(  I+K4  (D+E)  +K3A j 

1  +  R  (A+E)  +  K  m+E)  +  K,K  .E (A+D+E) 

3  4  3  4 


=  (MqV/aB  +  k2*FB)  P0c0  f  (2-44) 

max  2RT 

and 

AV3K3E  (l +K4  (A+D+E)  )  -  AV4K4E  (1+K3  (A+D+E)) 

AV  =  AV.  -  — - - - — - - - - 

1  1  +  K _  (A+E)  +  K  (D+E)  +  K -K .E (A+D+E) 

3  4  3  4 


Finally,  consider  the  coupling  between  the  two-step  and  exchange 
reactions.  The  system  will  be 


3 

A  +  B  +  E^C  +  D  +  E 

1  ++  2  ++  4  .  (2-45) 

F  +  B  t  C  +  B 

The  rate  equations  for  this  system  can  be  written 


+  |k  3+B+k  ^  +k  1  +  (A+E)j  AA  +  ^  -k^D^AC  +  ^  k  1  +A-K3+A-k^~cj  AD  =  0  , 


dAC 

dtT 


and 


dAD 

dt 


(k2++k3+B)AA  +  (k2~+k2++k3_D)AC  +  (k3+A+k3_c)AD  =  0  '  (2-46) 

|  k^+D-k  2+B j AA  -  k3~D  AC  +  ^k^”+k^+(D+E) +k  j~C+k  3+A j  AD  =  0  . 


If  both  steps  1  and  4  are  diffusion-controlled,  there  is  only  one  slow 
relaxation  and  the  parameters  are  given  by 
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k2+(R1  (A+E)+K1K4E(A+D+E))+k3+  (  K4DA+B ( 1 +K4 (D+E) ) +k 3_CK4D 

wr=k2+k3D  f+K,(A+E)  +  K,  (D+eT+~K,K  .E(A+D+E) 

7  4  14 

and  (2-47) 

AVt  (K1  (A+E)+K1K4E(A+D+E)J  -  AV4  (l+K!  (A+E) +K4E+K1K4E(A+D+E)) 

AV  -  AV3  l"+K  7  ( A+E )  +K  ,  (D+E )  +  K ,  K  E  ( A+D+E) 

1  4  1  4 

In  comparing  calculations  and  measurements,  frequency  (f)  will  be 
expressed  in  Hz,  where  f  =  oj/2h  .  Absorption  units  will  be  converted  from  Np 
to  dB  by  multiplying  by  8.7. 

2.3  RESONATOR  METHOD 

The  resonator  apparatus  is  illustrated  schematically  in  figure  2-4. 
The  resonator  used  in  these  experiments  is  a  72-liter  spherical  boiling  flask 
similar  to  the  others.  To  minimize  residual  losses,  it  is  suspended  in  a 
vacuum  chamber  by  loops  of  piano  wire  as  in  the  previous  experiments. 

The  flask  is  also  modified  with  a  small  ball-joint  vacuum  socket 
replacing  the  standard  neck.  The  mating  ball-joint  is  fitted  with  a  pair  of 
glass  tubes  which  allow  circulation  of  the  liquid.  Extending  the  input  tube 
a  few  centimeters  into  the  flask  produces  very  efficient  mixing  by  jet  action. 
The  exit  tube  is  flush  with  the  top  inside  surface  so  that  air  is  not  trapped. 
Flexible  tubing  is  used  to  connect  the  glass  tubes  to  a  corresponding  pair  of 
metal  tubes  which  pass  through  the  cover  plate  of  the  vacuum  chamber.  Addi¬ 
tional  flexible  tubing  is  also  used  to  connect  the  external  pump  and  reservoir 
to  the  system.  In  addition  to  storage  of  excess  liquid,  the  reservoir  serves 
for  adding  chemicals  which  are  then  circulated.  In  this  way,  the  surrounding 
vacuum  is  not  sacrificed  and  the  resonator  remains  physically  undisturbed. 

The  acidity  of  the  solution  is  continuously  monitored  by  means  of 
a  pH  meter.  Adjustment  to  a  desired  value  is  accomplished  by  adding  concen¬ 
trated  hydrochloric  acid  or  sodium  hydroxide.  Timing  tests  show  that  complete 
mixing  is  achieved  in  only  3  min.  with  the  pump  operating  at  2  liters/min. 
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2-4.  Schematic  of  Resonator  Apparatus 
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A  PZT  cylinder,  3.8  cm  in  length  and  diameter,  serves  as  a  trans¬ 
ducer  for  both  the  acoustic  measurements  and  degassing  the  liquid.  With  the 
flask  almost  completely  filled  with  distilled  water,  excess  dissolved  gas  is 
removed  by  boiling  at  reduced  pressure.  Cavitation  induced  by  driving  the 
cylinder  at  its  60  kHz  resonance,  helps  speed  up  the  process. 

The  reservoir  is  also  degassed  initially  at  reduced  pressure,  and 
also  while  dissolving  chemicals  which  may  contain  occluded  gases.  Accidental 
removal  of  desired  CO2  is  avoided  by  omitting  this  step,  unless  the  pH  is 
well  above  7  where  ionization  is  nearly  complete. 

In  the  measurement  mode,  the  liquid  is  at  atmospheric  pressure 
while  pressure  in  the  vacuum  chamber  is  reduced  to  less  than  1  mm  Hg.  The 
resonator  is  excited  at  a  selected  frequency  by  a  pulse  of  a  few  seconds  dur¬ 
ation  which  allows  the  sound  field  to  reach  equilibrium.  The  driver  used  for 
this  purpose  is  a  60  W  amplifier  capable  of  producing  100  V  rms  across  the 
load  up  to  frequencies  in  excess  of  100  kHz.  A  signal  generator  of  the  fre¬ 
quency  synthesizer  type  is  required  for  stability.  A  tone-burst  generator  is 
used  to  control  the  duration  of  the  pulse  and  the  time  delay  between  pulses. 

After  the  excitation  pulse,  the  transducer  is  automatically 
switched  from  transmit  to  receive  by  means  of  a  T/R  relay.  The  receiver 
system  consists  of  a  low-noise  preamplifier  and  a  tuneable  bandpass  filter. 
The  output  of  the  filter  is  fed  to  a  storage  oscilloscope  so  that  the  signal 
decay  can  be  monitored  during  the  tuning  procedure.  Measurement  of  the  decay 
rate  is  by  means  of  a  logarithmic  graphic  recorder  with  50  dB  range.  For 
convenience  in  comparing  resonator  measurements  with  sea  data,  the  measured 
values  are  converted  from  dB/sec  to  dB/km  by  dividing  by  the  nominal  sound 
speed  1.5  km/sec. 

Finding  the  resonances  is  a  very  tedious  process.  Because  the 
“sphere"  is  far  from  ideal,  calculation  of  resonant  frequencies  is  generally 
impractical  and  they  must  be  located  by  searching.  At  the  lower  frequencies, 
the  Q  may  be  so  high  that  the  frequency  cannot  be  stepped  more  than  1  Hz  at 


« 
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a  time;  otherwise,  a  resonance  may  be  missed  entirely.  At  higher  frequencies, 
bigger  steps  are  possible,  but  the  range  to  be  covered  is  greater  and  takes  a 
long  time  to  cover. 

Care  must  also  be  taken  when  adding  chemicals  because  the  sound 
speed  increases  causing  the  resonant  frequencies  to  increase.  Measurement 
of  the  shift  at  the  lowest  frequency  and  calculation  of  the  larger  shifts 
at  higher  frequencies  is  an  effective  means  of  following  the  changes. 

Certain  resonances  may  become  unusable  or  even  lost  and  others 
must  be  found  to  replace  them.  The  main  reason  for  this  difficulty  is  mode¬ 
coupling,  a  problem  that  has  not  been  extensively  discussed  in  the  earlier 
work.  When  resonances  involve  more  than  one  mode,  interference  effects  can 
be  troublesome.  If  the  frequency  spacing  is  great  enough,  the  decay  curve 
may  oscillate.  If  one  mode  then  decays  more  slowly  and  the  oscillations  die 
out,  the  data  may  be  usable.  However,  if  the  modes  are  closely  spaced  and  of 
comparable  decay  rate,  the  decay  will  appear  to  be  linear  but  the  slope  can 
change  with  small  changes  in  sound  speed.  Such  changes  in  residual  loss 
seriously  affect  the  accuracy  of  measurements,  particularly  at  the  lower 
frequencies  where  mode  interference  is  the  primary  source  of  error. 

Various  procedures  can  be  effective  in  alleviating  the  mode  prob¬ 
lem.  Absorption  can  be  increased  by  increasing  concentrations  up  to  the 
point  where  other  problems  arise  (precipitation,  polymerization,  etc.). 

For  pH-dependent  relaxations,  losses  may  also  be  increased  by  increasing  pH 
above  the  normal  range.  In  addition,  the  residual  losses  can  be  checked  by 
reducing  pH  to  the  point  where  the  chemical  losses  become  negligible.  If 
the  residual  loss  changes  excessively,  the  measurement  can  be  rejected  and 
the  procedure  continued,  altering  the  conditions  slightly  until  repeatable 
results  are  obtained. 

The  lowest  two  radial  modes  of  the  resonator  are  too  lossy  to  be 
of  use.  The  lowest  usable  resonance  in  all  cases  appears  to  correspond  to 
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the  third  mode.  For  the  72-liter  sphere,  the  corresponding  frequency  is 
nominally  12  KHz.  The  highest  usable  frequency  is  governed  by  viscous  losses 
which  increase  rapidly  above  100  kHz  making  the  resonances  difficult  to  find. 
The  same  problem  may  also  arise  at  the  lower  frequencies  if  the  chemical  loss 
becomes  too  large.  Combination  of  all  these  factors  results  in  a  measurement 
window  of  limited  range  in  both  frequency  and  absorption  coefficient.  The 
absorption  spectrum  must  fall  within  the  window,  if  the  relaxation  process 
is  to  be  analyzed. 

The  accuracy  of  the  resonator  method  is  difficult  to  access.  At 
the  lowest  resonance,  the  residual  loss  is  roughly  2  dB/km.  Incremental 
changes  as  low  as  0.2  dB/km  can  probably  be  measured  with  acceptable  accuracy, 
but  only  after  considerable  repetition  and  averaging.  Near  100  kHz,  residual 
plus  H2O  losses  increase  to  roughly  6  dB/km  and  the  minimum  limits  are  cor¬ 
respondingly  greater.  Estimates  of  relaxation  frequency  and  amax  involve 
averaging  over  the  10  -  100  kHz  range,  but  errors  at  the  higher  frequencies 
tend  to  dominate  and  the  accuracy  is  not  significantly  improved  over  the  low- 
frequency  value.  A  rough  estimate  of  the  rms  error  might  be  +_  10  percent 
for  amax  in  the  range  2-20  dB/km  and  also  for  fr  in  the  range  20  -  50  kHz, 
with  rapidly  increasing  error  outside  those  limits.  However,  possible  syste¬ 
matic  errors  in  the  procedure  may  tend  to  bias  the  results.  In  any  case,  the 
B(OH)3  values  for  nominal  sea  water  conditions  fall  below  the  limits.  There¬ 
fore,  measurements  are  made  at  increased  concentrations  and  compared  to  the 
sea  values  by  extrapolation. 


2.4  BORIC  ACID 

Simmons'  model  for  the  boric  relaxation  in  both  synthetic  sea  water 
and  NaCl  solution  is  the  two-step  ionization  equilibrium 


B (OH ) 3+OH 


(2-48) 


2-21 


TR-335N 


where  the  slower  step  2  controls  the  rate.  Argument  is  based  on  T-jump  meas¬ 
urements  of  relaxation  frequency  vs.  concentration  at  pH  =  8.0  (figure  2-5). 
Note  that  the  relaxation  frequency  is  always  below  the  resonator  range. 

If  ion-pairing  is  neglected,  the  data  of  figure  2-5  can  be  fitted 
by  equation  (2-36) 

k2+K1  [B(OH)3 
ur=  k2  +  1 +K 1  |B (OH) 3 

The  hydroxyl  term  [OH  ]  ~10-6  M  is  neglected.  Values  of  k2+  and  Ki  can 
be  determined  from  the  data  but  k2  is  too  small  to  measure  and  its  value 
must  be  inferred  from  the  overall  equilibrium  constant  Ki(1+K2)»  where 
k2  =  k2+/k2  • 

If  the  effects  of  ion-pairing  are  included,  the  equilibrium  of 
equation  (2-39)  shows  that  the  individual  association  constants  of  the  two 
borate  species  must  be  known. 

A  second  problem  is  what  values  to  use  for  calculating  AV.  This 
involves  not  only  the  association  constants  and  their  corresponding  AV's, 
but  also  pH.  The  reason  is  that  there  is  coupling  to  the  H2O  equilibrium. 

The  volume  change  for  the  H2O  equilibrium  at  25*C  is 


(2-49) 


AVw  =  V°(h+)  +  V°(0H~)  -  v(h2o|  =  -4.7  +  0.5  -  18.0  =  -22.2/kM  .  (2-50) 


Measurements  of  the  volume  change  for  the  boric  acid  equilibrium 
can  be  made  by  the  pH  vs.  pressure  method.  The  pertinent  equilibrium  con¬ 
stant  in  this  case  is 


[B(°H>4~]  H+ 

9 

(2-51) 

jBtOHTJ 
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where  [b  (OH)  4“j  =  |(b(OH)4")i  ]+  |  (b  (OH)  The  apparent  volume  change  is 

calculated  from  the  relation 


* 

AV 


-31nKB 

3P 


RT 


(2-52) 


The  equilibrium  constant  for  equation  (2-48)  is  K1(1+K2).  Thus, 

Kb  =  KWK i (1+K2) ,  where  Kw  =  [H+]  (OH”)  =  10~14,  and  equation  (2-52)  yields 


AV* 


AV  +  AV,  + 
W  1 


V^2. 

1  +K_ 


(2-53) 


The  estimated  values  for  aqueous  solution  at  25*C  are  AV*  =  -35.5/kM  and 
AVW  =  -22.2/kM.  With  K2  >>  1,  we  then  have  AV^  +  AV2  =  -13.3/kM. 

The  same  results  are  obtained  by  the  picnometer  method.  The  molal 
volumes  at  25*C  are  V°{B(OH)3)  =  39.2/kM  and  V°(b(OH)4  )  =  26.4/kM.  Given 
V°(OH~)  *  0.5/kM,  we  then  have 

AV1+AV2  =  V°(B(OH)4-)  -  V°(b(OH)3)  -  V°(0H-)  -  -13.3/kM  .  (2-54) 

Equation  (2-38)  with  K3  =  [H20]/lH+]  [OH-]  =  5.6  x  10* S/M  yields 


AV  i  | 

oh” 

)  + 

H+ 

AV  ~  AV  + 

2 

lH  + 

+ 

OH 

1 

l1  +  K1 

B  (OH)  3] 

1) 

Now,  for  pH  =  8.0  or  greater  [H+]  «  [OH_] ,  it  is  evident  that  coupling 
effects  to  the  H20  equilibrium  will  be  negligible.  At  pH  =  6.0  or  less 
AV  *  AV2  +  AVW,  but  absorption  is  negligible  in  that  range  anyway. 
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The  effects  of  NaB(OH)4  ion-pairing  on  absorption  also  appear  to  be 
negligible,  since  is  independent  of  NaCl  concentration  up  to  0.4M.  This 

shows  that  the  slow  step  is  unaffected  by  Na+  ion-pairing.  The  appropriate 
equilibrium  in  this  case  is  equation  (2-39)  with  the  values  K3  =  K5  =  2.1/M 
and  AV3  =  4V5  =  8/kM.  Then  AV  is  unaffected  despite  the  large  change  in  molal 
volume  for  NaB(OH)4°  formation. 

The  relaxation  frequency  obtained  from  equation  (2-40)  with 
[Na+]  >>  |b(OH)4"|,  is 

-  >2*  K,  [b|OH)3|  (l  *  KnaB|na*|)  .  ,2-561 

“C  ‘  k2  +  '  +  Kj  |B(OH)3|  (1  .  KNaB|».r|) 

The  only  effect  is  that  K]  p  +  KNaB  (Na+]j  appears  in  place  of  K-| .  This 
indicates  that  the  apparent  equilibrium  constant  of  boric  acid  in  NaCl  solu¬ 
tion  should  be  used  in  calculating  l<2  .  In  addition,  ion-pairing  of  NaOH° 
should  not  be  neglected  because  the  ionization  constant  of  water  changes 
significantly  in  NaCl  solution.  By  similar  arguments,  we  can  show  that  the 
apparent  equilibrium  constant  of  H2O  should  also  be  used  in  the  calculations. 
Table  2-2  lists  the  calculated  values  based  on  these  assumptions. 


Table  2-2.  B(OH)3  Rate  and  Equilibrium  Constants  Calculated  Using 

*  * 

Apparent  Constants  pKy  =  13.77  and  pKB  =  B.83  With  Estimated 
Diffusion-Control  Value  for  ki+  (Simmons  Theory) 


NaCl 

SSW 

K*  =  5.6  x  102/M 

K*  =  1.6  x  103/M 

+  10 , 

+  10 

k  ^  =3x10  /M  •  s 

k  =  3  x  10  /M*  s 

8 

7 

k  =  1.4  x  10  /s 

k  1  =  4.2  x  10  /s 

K2  =  154 

K2  =  53 

+  4 

k2  =  2.5  x  10  /s 

k2  =  2.2  x  10  /s 

2 

2 

k2  =  1.6  x  10  /s 

k2  =  4  x  10  /s 
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In  equation  (2-54),  AV'j  and  AV2  cannot  be  differentiated  by  absorp¬ 
tion  measurements  because  at  the  concentration  where  amax  can  be  measured 
accurately,  the  effect  of  AV^  is  very  small.  Simmons  proposes  that  the  inter¬ 
mediate  state  is  a  contact  ion-pair  with  negligible  volume  change  and  that 
all  of  the  AV  arises  from  the  second  step.  Ibis  is  consistent  with  other 
evidence  that  a  change  from  trigonal  planar  to  tetragonal  symmetry  is  involved 
in  the  ionization  of  boric  acid29  and  this  may  be  the  rate-controlling  process. 

Measurements  of  amax  vs.  concentration  are  compared  with  theory  in 
figure  2-6.  The  constituent  used  in  this  case  is  ^28407 • 1 OH2O  (borax).  The 
theoretical  curve  is  for  Simmons'  model. 

Figure  2-7  compares  measurements  of  amax  as  a  function  of  pH  for 
borax  and  boric  acid.  A  possible  reason  for  the  anomalous  pH  behavior  of 
borax  is  that  a  stable  aquo-complex  may  persist  which  inhibits  absorption  at 
high  pH  (equation  2-8).  The  pH  dependence  of  boric  acid  is  normal,  but  the 
magnitude  is  too  large. 

Despite  the  obvious  disagreement  between  theory  and  experiment, 
some  conclusions  can  be  made.  Note  that  the  measured  value  of  at 

pH  =  8.0  is  roughly  twice  that  of  sea  water  while  the  concentration  is  some 
20  times  greater.  At  lower  concentrations,  absorption  becomes  too  small  to 
measure,  except  at  higher  pH  values  where  the  lowest  measureable  concentration 
is  roughly  5  times  sea  water.  Extr apolat ion  to  sea  water  conditions  shows 
that  ijnax  is  more  than  an  order  of  magnitude  too  low. 


2.5  SYNTHETIC  SEA  WATER 

For  B(OH>3  absorption  measurements  in  synthetic  sea  water,  it  is 
desirable  to  eliminate  all  other  relaxations  that  cause  additional  absorption 
in  the  measurement  range.  This  can  be  done  by  o.  iting  certain  constituents. 
However,  coupling  between  equilibria  may  change  and  this  could  affect  the 
B(OH)3  relaxation. 
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Figure  2-6.  Resonator  Measurements  of 
in  Aqueous  B(OH)^  Compare' 
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Sulfate  will  be  omitted  with  no  apparent  effect  other  than  elim¬ 
ination  of  the  MgSO^  relaxation  which  would  otherwise  dominate  absorption 
in  the  10  -  100  kHz  range.  Magnesium  is  also  involved  in  the  MgCOj0  and 
MgB(OH>4+  relaxations,  both  of  which  produce  considerable  absorption  in 
this  range.  The  effects  of  omission  also  appear  to  be  small.  The  only 
constituents  in  addition  to  B(0H)3  that  are  necessary  to  approximate  the 
essential  sea  water  conditions  are  then  CaCl2  and  NaHCOj.  However,  NaCl 
may  be  useful  to  stabilize  carbonate  activity. 

The  T-jump  measurements  showed  that  NaHCO^  causes  an  increase  in 
the  B(0H>3  relaxation  frequency,  while  other  constituents  tend  to  reverse 
this  effect.  As  in  aqueous  solution,  the  relaxation  frequency  never  exceeds 
4  kHz  under  any  conditions. 

The  effects  of  CO2  on  the  B(OH>3  absorption  spectrum  in  0.4  M  NaCl 
solution  are  shown  in  figure  2-8.  At  the  highest  concentration,  the  measured 
relaxation  frequency  is  roughly  20  kHz.  Initially,  this  was  believed  to  be 
another  relaxation  because  the  relaxation  frequency  exceeds  the  T-jump  max¬ 
imum.  The  first  experiments  involved  NaCl  solution,  and  an  ion-pairing 
relaxation  was  suspected.  However,  attempts  to  match  the  spectra  with  a 
two-relaxat ional  model  were  not  successful.  Subsequent  experiments  showed 
high  relaxation  frequencies  without  additional  NaCl.  Here  the  results  are 
difficult  to  analyze  because  of  changes  in  carbonate  activity.  However,  the 
spectra  appear  to  be  almost  identical  in  both  0.1  M  and  0.4  M  NaCl  solutions. 
In  addition,  similar  results  are  observed  when  CaCl2  is  substituted  for 
NaCl.  The  observed  relaxation  is  clearly  the  sea  water  one,  the  higher  relax- 
at  ion  frequency  being  simply  the  result  of  increased  concentrations.  The 
discrepant  with  T-jump  measurements  can  be  explained,  if  pH  and  absorption 
are  controlled  by  different  relaxations  under  certain  conditions. 
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Figure  2-9  shows  the  absorption  spectra  with  CaCl2  substituted 
for  NaCl.  In  this  experiment,  the  NaHC03  concentration  is  reduced  to  0.4X 
(0.4  times  the  sea  water  value)  to  prevent  CaC03  precipitation  at  high  pH. 
The  absorption  spectra  are  then  measured  for  various  total  B(OH)3  concen¬ 
trations. 


In  both  experiments,  the  relaxation  frequency  is  approximately  a 
linear  function  of  the  varied  concentration.  Extrapolated  to  sea  water  con¬ 
centrations,  the  relaxation  frequency  is  approximately  equal  to  the  ocean 
value  f r ^  1  kHz. 

Now  amax  also  increases  linearly  with  either  concentration.  Values 
of  amax  in  NaCl  solution  corrected  for  concentration  and  pH  to  sea  water 
conditions  are  an  order  of  magnitude  too  small.  However,  with  CaCl2»  agree¬ 
ment  is  good.  It  is  apparent,  therefore,  that  the  omitted  constituents  play 
only  a  minor  role  in  the  B(0H)3  relaxation. 

Linear  dependence  of  relaxation  frequency  on  the  sum  of  concentra¬ 
tions  and  product  dependence  of  amax  indicates  the  exchange  equilibrium, 

B(OH)3  +  co32_  +  h2°  t  B(OH)4“  +  HC03'  .  (2-57) 


For  the  exchange  equilibrium,  we  have  from  equation  (2-42) 


“r  =  k+  ( 


2- 

i\  -/ 

-  1  f 

3 (OH) 3 

+ 

C°3 

1)  *k  ( 

B  ( OH )  4  +  HC03 

(2-58) 


max 


(AV) 2k+  |b(OH) 3|  |C032")  p0C0 

-  — -  ' 
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The  reason  for  the  negligible  effect  of  Na+  ion-pairing  is  appar¬ 
ently  the  same  as  proposed  for  aqueous  solution,  i.e.,  the  rates,  Av's  and 
association  constants  are  similar  for  both  reaction  paths.  The  suggested 
exchange  scheme  in  NaCl  solution  is 


Na+  B(OH)3  +  CO  2"  +  H20 

3  +t 


Na  +  B (OH) .  +  HCO- 

4  3 

4  ++ 


(2-59) 


B (OH)  +  NaCO,  +  HO  t  NaB (OH)  +  HCO. 
3  3  z  4  3 


For  [Na+]  » 


B (OH) 4 


C03z_  |  in  equation  (2-44), 


.  +/  B(OH)3] 

id  I - —  + 

r  1  *  1  +  K3  ]Na+ 


CO. 


B  (OH) 


hc03~ 

1  \ 

1 

f  K4 

Na+ 

+  k. 


/K3_  1 

Na+| 

|  B (OH)  3j 

l  ' 

+  K3 

Na+ 

+  CaC03  |  + 


K3  Na 


B  (OH) 


HCO- 


(2-60) 


1  +  K4  |Na 


av  =  av 1  - 


AV3K3  Na 


1  +  K3  |Na 


AV4K4  ! 

Na+| 

1  +  K4 

Na 

1 

In  this  case,  the  forward  rate  is  independent  of  [Na+] ,  if 
ki+  =  k 2+ -  The  reverse  rates  cannot  then  be  equal  because  K3  >  K4,  and 

k)  /k2  3  K4/K3.  For  calculations,  the  values  used  are  Ki  *  4,  K3  =  1.6/M, 

*  *  ,  , 

K4  =  2.1/M,  pKB  =»  9.0,  and  pK2c  =  9.6.  Assuming  k^  =  k2  ,  the  relaxation 

frequency  can  be  written  u>r  *  k+  X  effective  concentration.  The  value  of  k+ 

estimated  by  fitting  the  data  of  figure  2-10  is  7.5  x  106/M*s. 
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In  figure  2-11,  the  data  are  plotted  vs.  the  concentration 

product  |b(0H)3|  x  j Z C03^—  j .  The  value  of  Av  is  calculated  by  curve-fit. 

In  NaCl  solution,  Av  =  -7/kM.  The  value  for  the  exchange  is 

AV  *  VO(b(0H)4~|+VO(hC03")-VO(b{0H)3)-VO(C032")-v(H20)  =  -8.2/kM.  (2-61) 

In  equation  (2-60),  we  then  have  AVi  =  -8.2/kM,  AV3  =  12/kM,  and  AV4  =  8/kM. 

The  calculated  value  is  then  Av  =  -7.4/kM,  which  is  in  good  agreement  with 
the  measured  value. 

In  CaCl2  solution,  the  curve-fit  value  is  Av  =  -21/kM.  The  scheme 
for  CaCl2  is 


Ca2+  +  B (OH)  +  C0,2~  +  H_0  1  Ca2+  +  B(OH)  "  +  HC0,~ 

3  3  2  4  3 

3  ++  2  4  +*  •  (2-62) 

CaCO.°  +  B (OH)  +  HO  t  CaB(OH) .+  +  HCO  ~ 

3  3  2  4  3 

In  this  case,  the  approximation  jca2+]  »  Jb(OH)4  ]»  [co32-]  *s  not  accurate 

and  the  exact  equations  are  used.  For  calculation  purposes,  the  values 

*  * 

are  pKB  =9.1,  pK2c  =  9.7,  K4  =  26/M  and  K3  =  125/M.  For  Av4  =  0  and 
AV3  =  12.5/kM,  the  calculated  value  is  AV  =  -16/kM,  which  is  significantly 
too  low  in  magnitude.  The  likely  explanation  is  that  K3  is  much  larger  than 

the  sea  value  at  the  lower  ionic  strength,  making  AV  =  AV-j  -  AV3  =  -20. 7/kM, 

which  is  in  good  agreement  with  the  measured  value.  In  any  case,  the  strong 
effect  of  Ca2+  is  clearly  due  to  the  asymmetry  in  the  association  constants 
and  Av' s  of  the  two  Ca2  +  ion-pairs. 

In  figure  2-12,  the  pH  dependence  of  is  compared  with  theory. 

From  equation  (2-58),  we  have 


+ 

a  *k 
max 


[B] 


1_ 
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where  [B]  and  ICO2)  are  total  concentrations,  x  =  10PH-PKB,  and  y  =  ioPh-Pk2C. 
Agreement  is  taken  as  further  support  of  the  exchange  equilibrium  as  the 
dominant  mechanism  in  sea  water. 


In  sea  water,  the  effects  of  Na+  tend  to  cancel  out.  The  same 
appears  to  be  true  of  Mg2+.  However,  it  is  more  difficult  to  demonstrate 
because  two  additional  ion-pair  relaxations  must  be  considered.  Figure  2-13 
shows  the  absorption  spectrum  of  the  B (OH) 3/C02/Mg/Ca  system  at  pH  =  8.5. 

All  concentrations  are  sea  water  values,  except  B(0H)3,  which  is  4.6X. 
Reduced  NaHC03  concentration  is  no  longer  necessary  to  prevent  precipitation 
of  CaC03  because  competitive  MgC03°  ion-pairing  keeps  the  concentration  from 
becoming  excessive.  The  solid  lines  represent  the  estimated  spectra  of  the 
three  component  relaxations  B(0H)3,  MgC03°  and  MgB(OH),j+  (reference  23). 

Figure  2-14  shows  the  calculated  pH  dependence  of  the  three  com¬ 
ponents.  The  data  points  are  12  kHz  measurements  with  the  other  components 
subtracted  to  estimate  the  B(0H)3  values.  The  curve  fit,  extrapolated  to 
sea  water  concentration,  is  about  25  percent  too  high,  which  indicates  that 
Mg  causes  an  increase  in  amax.  A  detailed  analysis  of  the  interaction  is 
difficult  because  the  MgB(0H)4+  relaxations  are  also  slow  and  coupling 
effects  are  not  easily  estimated  when  the  rates  are  comparable. 


In  figures  2-15  and  2-16,  the  pH  dependence  of  the  full  SSW  system 
is  examined  at  two  widely  different  temperatures.  Analysis  of  sea  data 
(reference  18)  showed  that  B(OH)3  absorption  data  can  be  expressed  in  terms 
of  fr  as  a  function  of  temperature  only,  and  (al)max  as  a  function  of  pH 
only.  Thorp's  values  at  T  »  4*C  and  pH  =  8.0  are  used  as  reference  and  the 
empirical  formula  can  be  written 


f 

r 


1  x  10 


(T-4 ) /I °°kHz 


a 

max 


0.1  x  10  (pH_8)fr  dB/km 


(2-64) 
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where  the  nominal  pH  range  is  7.8  -  8.2.  Linear  regression  analysis  of  the 
same  data  by  Schulkin  and  Marsh313  and  by  Lovett31  yielded  a  coefficient  for 
the  pH  exponent  somewhat  less  than  unity  which  would  account  for  small  devi¬ 
ations  from  the  asymptotic  value  in  this  range. 

For  comparison  of  the  formula  values  with  the  results  of 
figures  2-16  and  2-17,  the  coefficients  are  multiplied  by  4.6  to  correct  for 
the  increased  B(0H)3  concentration.  At  9 °C  and  pH  =  8.0,  the  curve-fit  value 
of  figure  2-16  is  0.57  dB/km  and  the  model  value  is  0.52  dB/km.  At  27*C,  the 
corresponding  values  are  0.78  and  0.79  dB/km.  Obviously,  such  excellent 
agreement  is  somewhat  fortuitous,  since  errors  of  more  than  10  percent  might 
easily  be  expected. 

The  linear  approximation  for  (a\)max  vs.  1 0PH  becomes  inaccurate  at 
high  values  of  pH.  In  figure  2-17,  the  sea  data  and  resonator  data  are  com¬ 
pared  to  the  actual  theoretical  curve.  The  inflection  point,  depends  on  the 

*  ★ 

values  of  pKg  and  pK2c*  both  of  which  depend  on  temperature.  However,  for 
calculation  purposes,  25 ‘C  values  are  appropriate  for  the  resonator  measure¬ 
ments,  while  errors  in  the  sea  water  temperature  range  are  evidently  not 
significant  because  AV  changes  evidently  cancel  out  effects  of  pK  changes. 
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Section  3 

DISCUSSION 


The  exchange  model  appears  to  give  good  account  of  B(0H>3  absorp¬ 
tion  in  synthetic  sea  water  and  the  ocean  measurements  as  well.  However, 
there  are  still  a  numbe'  of  discrepancies  relating  to  the  T-jump  measure¬ 
ments.  Note  that  the  exchange  of  a  hydroxyl  ion  does  not  involve  pH  change. 
Thus,  it  is  only  through  other  coupled  equilibria  that  absorption  and  T-jump 
measurements  are  related. 

In  figure  3-1,  the  relaxation  vs.  concentration  predictions  of 
exchange  theory  are  compared  with  Simmons'  SSW  results.  For  the  exchange 
case,  relaxation  frequency  is  a  linear  function  of  BfOHJj  concentration. 

The  finite  value,  when  the  concentration  approaches  zero,  is  due  to  HCO3-, 
which  is  the  predominant  CO2  species  at  pH  =  8.0.  T-jump  data  agree  within 
experimental  error,  except  for  the  highest  value  where  it  appears  to  approach 
a  limit  as  in  two-step  theory. 

In  figure  3-2,  the  relaxation  frequency  vs.  pH  measurements  in  sea 
water  by  Yeager,  et  al  are  compared  with  the  two  theories.  At  low  pH,  the 
relaxation  frequency  approaches  a  constant  in  either  theory.  At  high  pH,  the 
relaxation  frequency  for  exchange  theory  increases  approaching  another  con¬ 
stant  value  determined  by  CC^-  species.  Two-step  theory  predicts  a  decrease 
as  the  B ( OH ) 3  concentration  decreases.  The  data  are  in  reasonable  agreement 
with  either  theory,  except  at  the  highest  value,  where  the  two-step  seems  to 
give  the  correct  trend. 
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an 


Relaxation  Frequency  (kHz) 


(B  (OH)  3]  (mM) 


Figure  3-1.  Simmons'  T-jump  Measurements  of  B(OH)^ 
Relaxation  Frequency  vs.  Concentration  in  SSW 
Compared  With  Two-step  and  Exchange  Theories 
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Figure  3-2.  T-jump  Measurements  of  B(OH)3  Relaxation 
Frequency  vs.  pH  in  Sea  Water  Compared  With 
Two-step  and  Exchange  Theories 
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A  coupled  system  for  the  two  reactions  in  sea  water  would  be 


h2o  +  b(oh)3  +  oh"  +  co32 
t 

h2o  +  (b(0H)4")1  +  co32' 


3 

2 


+  HC03  +  OH 
4  it 

+  co32"  +  H20 


(3-1) 


where  step  3  is  the  exchange  and  1-2  is  the  two-step. 


In  a  four  state  system,  there  are  three  eigenf requencies.  If 
steps  1  and  4  are  both  fast  compared  with  steps  2  and  3,  two  are  fast  and 
the  other  is  slow.  Therefore,  pH  and  absorption  would  follow  the  same  slow 
relaxation,  regardless  of  whether  step  2  or  3  controls  that  rate.  The  "short 
circuit"  would  be  eliminated  if  either  steps  1  or  4  were  slow  also.  There 
is  no  evidence  for  a  slow  relaxation  in  the  carbonate  equilibrium  (step  4). 
However,  step  1  can  be  slow,  if  the  0”  association,  like  the  exchange,  is 
not  diffusion-controlled  as,  Simmons  assumed. 


Since  k2+  places  an  upper  limit  on  the  two-step  relaxation  fre¬ 
quency,  it  is  clear  that  the  new  value  must  ce  approximately  the  same  as 
before.  However,  to  fit  the  boric  acid  absorption  data  of  figure  2-7  with 
the  values  AV^  *  0  and  AV2  =  -13.3/kM,  a  much  higher  reverse  rate  is 
required,  i.e.,  k2~  ■  7.5  x  10^/s.  The  other  values  must  be  determined  by 
curve  fit  to  T-jump  data.  Approximate  methods  are  no  longer  valid  and  the 
eigenf requencies  must  be  calculated  from  the  roots  of  the  determinant  of 
equation  (2-35) 


+ 


(3-2) 


Curve  fit  of  the  lower  mode  II  to  the  NaCl  T-jump  data  of  figure  2-5  is  shown 
in  figure  3-3.  The  values  of  the  resulting  rate  and  equilibrium  constants 


Figure 


-3.  T-jump  Measurements  of  B(OH)3  Relaxation  Frequency 
in  NaCl  Solution  Compared  With  Present  Theory 
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ace  given  in  table  3-1.  Note  that,  although  the  higher  mode  I  falls  in  the 
experimental  range,  the  effects  may  not  be  noticeable  if  the  excitation  is 
small. 

Table  3-1.  B(0H)3  Rate  and  Equilibrium  Constants 

in  NaCl  Solution  (Present  Theory) 


In  the  normal  mode  model  for  aqueous  BJOHJj,  the  two  asymptotic 
equilibria  ace 

B (OH)  0H_  2  ^B(OH)  ~  J  x 


and 
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However,  the  relaxation  frequency  w1  is  not  always  greater  than  <*>2  as  it 
would  be  in  the  diffusion-controlled  case.  For  the  coupled  system,  the  two 
conditions  are  given  by  the  following  approximation: 


taJ 


0) 


II 

I  : 


B (OH) . 


k1  K2 
1+K_ 


<*>1  « 


k2+K1  B(0H)J 
“il'*’  k2  +  1  +  K1  |B  (OH)  3| 

+  K1  +  |b(0H)3  +  k2“ 


(3-4) 


Simmons'  measurements  show  no  measureable  changes  with  temperature  which 
indicates  that  all  temperature  variation  in  the  equilibrium  constant  comes 
from  k2~.  This  conclusion  is  also  supported  by  the  rather  large  temperature 
dependence  of  o,,^.  Excitation  of  step  1,  via  step  2  is,  therefore,  small 
for  >>  o»2  and  the  lower  mode  dominates  pH  change  throughout. 

In  the  case  of  amax,  the  upper  mode  is  directly  excited  through  AV2 
at  small  concentrations.  The  resulting  absorption  is  small  and  the  relaxa¬ 
tion  frequency  is  still  well  below  the  range  of  present  equipment.  At  higher 
concentration,  the  upper  mode  is  not  excited  because  AV^  *  0  and  the  lower 
mode  dominates. 

In  the  SSW  model,  only  step  4  is  assumed  to  be  diffusion-controlled. 
Then  the  rate  equation  (2-45)  can  be  approximated  to  give  the  following 
determinant: 
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Thus,  there  is  a  rapid  and  complete  mode  change.  In  this  case,  pH  change 
follows  the  lower  mode,  while  the  absorption  dominated  by  the  exchange  path 
branches  to  the  upper  mode.  The  short-circuit  condition  mentioned  earlier  is 
clearly  governed  by  k^  .  For  the  diffusion-controlled  case,  kj  >>  k2  +  k2+* 
and  the  transition  would  be  above  the  range  of  interest. 

The  present  model  does  not  resolve  the  pH  discrepancy  in  figure  3-2, 
since  the  relaxation  should  follow  the  exchange.  However,  at  pH  =  9.5,  sea 
water  is  highly  supersaturated  with  CaC03  and  precipitation  could  account  for 
the  reduced  value. 


Figure  3-4.  T-jump  Measureme 
Frequency  in  SSW  Compared 
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Section  4 

CONCLUSIONS 

— ^  Resonator  measurements  indicate  that  the  principal  BfOH)^  absorp¬ 
tion  mechanism  in  sea  water  is  an  exchange  between  the  B (OH) 3/B (OH) 4  and 
HC03~/C032-  equilibria  in  which  coupling  to  the  Ca^-  ion-pairing  equilibria 
provides  most  of  the  molal  volume  change.  If  either  Ca  or  CO2  are  omitted, 
B(0H>3  absorption  is  an  order  of  magnitude  too  low. 

T-jump  measurements  indicate  a  two-step  B(0H)3  ionization  process 
in  H2O.  In  sea  water,  both  relaxations  occur  and  are  coupled.  Present  theory 
shows  that  pH  and  absorption  may  have  different  relaxation  frequencies  as 
experimental  data  indicates. 

The  empirical  pH-dependent  absorption  formula  derived  from  ocean 
measurements  gives  predictions  that  are  consistent  with  the  exchange  relaxa¬ 
tion  model. 
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